Effect of initial correlations on short— time decoherence 
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We study the effect of initial correlations on the short-time decoherence of a particle linearly 
coupled to a bath of harmonic oscillators. We analytically evaluate the attenuation coefficient of a 
Schrodinger cat state both for a free and a harmonically bound particle, with and without initial 
thermal correlations between the particle and the bath. While short-time decoherence appears to 
be independent of the system in the absence of initial correlations, we find on the contrary that, 
for initial thermal correlations, decoherence becomes system dependent even for times much shorter 
than the characteristic time of the system. The temperature behavior of this system dependence is 
discussed. 
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Environment induced decoherence plays a fundamental 
role in many areas, ranging from quantum cosmology 0] 
and the theory of quantum measurement j2| to quantum 
information and quantum computing Experimental 
investigations of the decoherence process have recently 
been reported in Refs. Q-^. Environmental decoherence 
can be defined as "... the (irreversible) loss of quantum 
coherence of a quantum system due to its coupling to 
an environment" Q ; it manifests itself in the dynamical 
suppression of interference phenomena [^,^. A simple 
way to quantify the destruction of coherence is thus to 
consider the superposition of two localized wave pack- 
ets (Schrodinger cat state) and to look at the decay of 
the interference term, as measured for instance by the 
attenuation coefficient [see Eq. (|ll]) below]. An impor- 
tant characteristic of decoherence is that it occurs on a 
very fast time scale, usually much shorter than the en- 
ergy dissipation scale. We note that the short-time limit 
of decoherence has recently attracted a renewed inter- 
est in the literature [p^| - |l4| . Interestingly, Braun, Haake 
and Strunz |^,^ have identified a new regime of fast 
decoherence beyond the usual Golden Rule regime. In 
this limit of large separations between the wave packets 
(interaction-dominated decoherence), quantum decoher- 
ence appears independent of both the system and the 
heat bath. 

Most studies of environment induced decoherence 
make use of the simplified assumption that the system 
and the bath are initially uncorrelated Then, the 

initial composite density operator for system plus bath 
can be factorized into a product of a system operator and 
a bath operator. However, in the general, and more re- 
alistic case, initial correlations can be present and the 
latter might affect the decoherence process, especially 
at very short times as first pointed out by Romero and 
Paz In fact, as recently shown by Ford, Lewis and 
O'Connell initial thermal conditions can dramat- 

ically modify the decoherence rate. They find that the 
decoherence time for an unbound particle becomes in- 
dependent of the strength of the dissipation when it is 



initially at the same temperature as the bath. Our aim 
in this paper is to complement the discussion of the ef- 
fect of initial thermal correlations on short-time decoher- 
ence by analytically comparing the attenuation factor of 
a free particle with that of a harmonic oscillator, with 
and without thermal initial correlations. As our main 
tool, we shall employ the by now standard model of a 
system linearly coupled to a bath of harmonic oscillators. 
This model is exactly solvable and the exact reduced den- 
sity operator of the quantum system is conveniently ob- 
tained within a path integral approach [ pT[ . This model 
has been extensively used in decoherence studies for the 
case of factorizable initial conditions, building up on the 
work of Caldeira and Leggett It has been extended 
later on to correlated thermal initial conditions by Hakim 
and Ambegaokar for a free particle and by Morals 
Smith and Caldeira for a harmonically bound particle 
|20,pT| (see also the work of Grabert and collaborators 
|22|). In the following we shall use the exact results 
derived by Morals Smith and Caldeira ^l[] to compute 
the time evolution of the superposition of two Gaussian 
wave packets separated by a distance 2a and of width 
b, ipix) = ci exp[-(a:; - a)'^/Ab'^] + c-i cxp[-(x + af/Ah'^]. 
The corresponding initial density operator, ps{x,y,0) = 
'ip{x)'ip*{y), is given by 



Ps{x,y,Q) = Icipexp 
+ |c2|^exp 
-I- C1C2 exp 
-I- c[c2 exp 
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We calculate the short-time limit of the attenuation co- 
efficient for a harmonically bound particle with and with- 
out initial thermal correlations and compare the obtained 
results with those of a free particle. We find that for ini- 
tial decorrelation, the attenuation coefficient is identical 
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for both the free particle and the hnear oscillator, that 
is to say, decoherence is independent of the nature of the 
system. On the other hand, for initial thermal correla- 
tions, we find that the attenuation coefficient acquires 
a system dependent part, even for times much smaller 
than the characteristic time of the system. This shows 
that in the presence of initial thermal correlations deco- 
herence is system dependent. At very high temperature, 
this system dependence turns out to be negligibly small. 
However, in the opposite limit of low temperature, it be- 
comes increasingly important. 

The model. We consider a particle of unit mass, mov- 
ing in an external potential V{x) and linearly coupled 
through its position to a set of independent harmonic os- 
cillators with mass rrii and frequency tOi. The composite 
Hamiltonian is written in the form 



H = 



Pi 
2 



V{x) + ^ X dxi 



E 



2m, 



1 



2 2 



f<2 
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where the C^'s are coupling constants. In the limit of 
infinitely many oscillators, the bath is entirely character- 
ized by the spectral density function, 

I{lo) = J y 5{lo - uJi) = 27W 9{w - cu,) , (3) 

i 

where the last equality defines Ohmic damping with 
damping coefScient 7. Here lUc is a cutoff frequency that 
will be replaced by infinity in all convergent integrals. 
The reduced density operator psit) of the system at time 
t is obtained after tracing out the bath degrees of free- 
dom. In coordinate representation it reads 

ps{x,y,t)^ J dx'dy'dRdll'dQ' K{x,R,t;x',R',0) 

xK*{y,K,t;y',Q',0){x,Il'\p{0)\y',Q') . (4) 

Here R, R' and Q' collectively denote the coordinates of 
the bath, K is the propagator and p(0) the initial den- 
sity operator of the composite system. If we assume that 
at t = the system and the bath are uncoupled and 
that the latter is in thermal equilibrium, then the initial 
density operator can be written as po{x' ,11' ;y' ,Q') = 
ps{x' ,y' ,0)peq{'R' ,Q'), where ps{x',y',0) describes the 
initial state of the system and peg(R',Q') is the equi- 
librium density operator of the bath. On the other 
hand, if the system and the bath are initially coupled 
and in thermal equilibrium with each other, p{0) can- 
not be factorized into a product of a system and a bath 
operator anymore. Instead, we have pq(x' ^^''^y' ,Q') = 
Ps{x',y',0)peq{x\y','R',Cl'), where peq{x' ,y' ,R' ,Q') is 
now the equilibrium density operator of the composite 
system. We shall refer to these two cases as (i) uncorre- 
lated initial conditions and (ii) thermal initial conditions, 



respectively. We mention that other correlated initial 
conditions have also been considered (see Refs. [p0|-p2t). 
The integrals over R, R' and Q' in Eq. (||) can now be 
performed exactly, yielding 

Psiq,C,t)^ fdq'deJ{qA,t;q\C,0)Ps{q',e,0), (5) 



where we have introduced the center of mass and rel- 
ative coordinates q — {x + y)/2 and ^ ^ x — y. The 
dynamics of the particle is completely determined by the 
propagating function J(g, ^, i; q', 0). Equation (^) can 
be considered as the full solution of the master equation 
describing the time evolution of the dissipative particle 
with (and without) thermal initial conditions. In the di- 
agonal case ^ = 0, the propagator J is given by (we put 
h — kB = l throughout the paper). 



J{q,0,t;q',e,0) 



m 

2tt 
X exp 



exp 



•,(^a{t)q'(' -N{t) q^' 



eq- 



with 



and 



a{t) = K{t) -h 7 - 2eE{t) , 
A(t) = C{t) - eE{tf , 

^ 2k ' 



(6) 
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The exact expressions for the coefficients K{t), N{t), 
C{t) and E{t) have been derived in Ref. ||2l[ for a linear 
oscillator with frequency ojq initially in thermal equilib- 
rium with a heat bath at inverse temperature /3 — 1/T. 
They are reproduced in the appendix for completeness. 
Remarkably, the form of the propagating function (|^) 
remains the same for uncorrelated initial conditions, as 
well as for a free particle. The case of uncorrelated initial 
conditions is recovered by setting e to zero and keeping 
only the first double integral in Cuj{t), Eq. (25), [K{t) 
and N{t) being unchanged], while the unbound parti- 
cle is obtained by letting the frequency wq go to zero 
The quantity k that appears in Eq. (9) is equal 
to the variance of the position {x'^)eq in equilibrium [see 



Eq. (|2|) in appendix A]. The factor e (= 1/2A in the 
notation of |2l|), which stems from the initial thermal 
correlations, will turn out to be important in the follow- 
ing discussion. Its asymptotic behavior at high (T^ujq) 
and low temperature {T <ti luq) is respectively given by 
£ ~ (/3wo)^/2 <C 1 and s ~ /Sujo > 1. The diagonal den- 
sity operator of the system at time t can then be easily 
obtained by combining Eqs. (0), (||) and (^. We find 



ps{q,0,t) - exp 
-I- exp 
-|- 2 exp 



2PiNq - aaf + 0^(1 + Sb^A)^ 



2h'^{Nq + aa)^ + 0^(1 + 8b^A)e 



a^{AA + 2b^a^) + 2b'^N^q^ 



2 



X COS 



2aNq{l + 2b^s) 



(10) 



where we have defined ct^ = (1 + 86^ A) (1 + 2b'^e) + ib'^a^ 
(for simplicity, we have put ci equal to C2). Equation(|Tc|) 
is written as a sum of three terms. The first two terms 
correspond to two separately propagating wave packets, 
while the third one, containing the cosine, is an interfer- 
ence term. The attenuation coefficient a{t) is defined as 
the ratio of the factor multiplying the cosine to twice the 
geometric mean of the first two terms. It follows from 
Eq. (O) that 



-4a2A + a2(l + 862A)e 



(11) 



a(t) = exp 



This expression is still exact. The attenuation factor(|l] 
is the measure of decoherence we shall use in what follows 
to investigate the short-time limit of decoherence. To be 
more specific, we shall place ourselves in the limit of weak 
coupling between the system and the bath, 7 <C wp, and 
assume that the time is small compared to the relaxation 
time, <^ 1 and also small compared to the evolution 
time of the harmonic oscillator, toot <^ 1. 

Uncorrelated initial conditions. We begin by con- 
sidering initial decorrelation between the system and the 
heat bath. Physically, this corresponds to the situa- 
tion where the system and the bath are isolated prior 
to their coupling att = 0, the system being effectively at 
zero temperature. In this case £ = 0, a{t) = N{t) and 
A{t) = C{t) (with only the first double integral). In the 
high-temperature limit, T3>7,Ci;o, we approximate the 
hyperbolic cotangent in Eq. ( p^ ) by coth x ~ Ex- 
panding a{t) and A(t), Eqs. (7) and (8), in lowest order 
in time, we find that the attenuation coefficient (11) is 
given by 



"0 



it) — exp 



for the free particle (i<C7 ^), and by 



ag"{t) = exp 



1264 ^ sb'^Toj^^t^ + 3^2 



(12) 



(13) 



for the harmonic oscillator (< ^ ^o^)- For very short 
times, we thus obtain the same attenuation coefhcient 
for both the free particle and the linear oscillator. 



ao{t) ~ exp 



364 



(14) 



Interference patterns between the two superposed wave 
packets are hence destroyed according to a stretched ex- 
ponential on a time scale, td = {3b^ /jTa'^Y^^ . The deco- 
herence time T£) depends solely on the friction strength 
7, the temperature T and the parameters of the initial 
wave packets (separation a and width 6) , and not on any 



system specific quantity, like the frequency of the oscil- 
lator for example. Since we look at times much shorter 
than the characteristic time of the system, T5 ^ ^0^, 
such a system independence of the decoherence time is to 
be expected. However, as we shall see, this is only true 
for uncorrelated initial conditions. This result is remi- 
niscent of the universal regime discussed by Haake and 
coworkers, where the decoherence rate was also found 
to be independent of the nature of the system p^ , p^ . 
However, the short-time attenuation factor (|l^) does not 
quite belong to this regime of very fast decoherence char- 
acterized by a Gaussian decay law, exp[— (i/r^i)^]. On 
the other hand, the cubic time dependence clearly indi- 
cates that expression (14) goes beyond the (long-time) 
Golden Rule regime, t^ujQ^, and its typical exponential 
decay, exp[— t/r£)]. Here we have a much slower initial 
decoherence compared to the Golden Rule expression. A 
similar short-time cubic dependence of the decoherence 
factor already appears in Refs. ||l8| and We also 

mention that the short-time approximation used in the 
derivation of Eq. (|lj) amounts to neglecting the spread- 
ing of the wave packet that appears in the denominator of 
the attenuation coefhcient. For instance, the short-time 



spreading of the free wave packet in Eq. (12) is given by 
(Aa;2(t)) = 6^ t'^/Ab'^ + 2Tjt^/3, which reduces to the 
initial width 6^ for very short times. 

Thermal initial conditions. We now turn to the 
case where the system and the bath are initially corre- 
lated and in thermal equilibrium. First, it should be 
realized that thermal initial conditions not only modify 
the coherence time of the system, they also directly affect 
the overall coherence length As a matter of fact, we 
easily see from the coordinate representation of the free 
particle equilibrium density matrix, {x\ exp[— /3p2/2] \y) ~ 
exp[— (a; — y)'^/2(3], that there is an exponential cutoff in 
x — y over distances of the order of the thermal de Broglie 
wavelength, A ~ 1/Vt. The two wave packets can there- 
fore only coherently interfere if they are separated by a 
distance smaller than A. This distance is extremely small 
at high temperature, but becomes appreciable when the 
temperature is very low. 

In the limit of high temperature, we can compute the 
attenuation factor ( |ll| ) for the free particle, for times 
smaller than the relaxation time, t ^ 7~^, by expanding 
the functions a(t) and A{t) up to lowest order in t, in 
analogy with the previous section (now keeping all the 
terms). For very short times, this leads to 

2a^Tt^ 



exp 



exp 



464 + ATbH"^ 



t2 



264 



(15) 



The short -time expression dTs] ) is equivalent to the result 
recently obtained by Ford, Lewis and O'Connell using an 
exact method based on quantum distribution functions 
pl|. It is worth noticing that Ford and O'Connell have 
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also derived Eq. ( |l5| ) with a more elementary method 
making only use of basic quantum mechanics and equi- 
librium statistical mechanics | [lO| . This approximate 
method is valid in the limit of vanishingly small fric- 
tion. Equation ( p^ shows that, in the presence of initial 
thermal correlations between the system and the bath, 
the short-time dependence of the exponent of the de- 
coherence coefficient is now quadratic. This has to be 
contrasted with the cubic dependence obtained for ini- 
tial decorrelation in Eq. (p^). Moreover, the decoherence 
time, td = (26*/Ta^)^/^, appears to be independent of 
the friction coefficient 7. This remarkable result indicates 
that at high temperature, decoherence can occur without 
dissipation 

Similarly, the attenuation coefficient for the harmonic 
oscillator can be calculated in the limit t ^ f^cT^- Still in 
the high-temperature limit, we find that for very short 
times it is given by 



exp 



exp 



(1 + 262£)1 , 



(16) 



with £ ~ (/3a;o)^/2. Contrary to Eq. (|T^), we observe that 
for initial thermal conditions, the short-time attenuation 
coefficients for the free particle (^5|) and the harmonically 
bound particle ( p^ are not identical, even for times much 
smaller than the characteristic system time rg . Equation 
( |l6| ) indeed contains an additional, system specific cor- 
rection, that depends on the frequency of the linear oscil- 
lator Wo, and on the temperature T. This term finds its 
origin in the initial thermal correlations existing between 
the system and the heat bath at t = (see the discussion 
below). In the limit of high temperature, T ^ ujq, this 
correction, which is proportional to e, is negligibly small. 
However, as we shall see next, it becomes increasingly 
important as the temperature is lowered. 

In the low-temperature limit, we replace the hyper- 
bolic cotangent in Eq. ( p^ by unity, coth a; ~ 1. For 
very short times, we find that the attenuation coefficient 
for the free particle reads (t<C7~^), 



<T(i) 



exp 



exp 



9 



(ln7t-|-7e -3 + ln2)<^ 



u 2 



(17) 



where 7e is Euler's constant. We note that in contrast 
to the high-temperature expression (15), Eq. (^7|) now 
explicitly contains the damping coefficient 7 (see the re- 
cent discussion in Ref. ||2^ on this point). We also note 
that the exponent in Eq. (^) is actually negative since 
7t < 1. The "i^lni" behavior of the decoherence factor 
in Eq. (^7|) is consistent with the result found by Romero 
and Paz for the superposition of two translations Eq] . In 



an analogous way, we find that the short-time expres- 
sion of the attenuation coefficient for the linear oscillator 
is given by {t<$:ujQ^), 



~ exp 



(lnwo< + 7e)^^(l + 26^e) 



20^7 
r2a2 

—jt^ lnujotil + 2b^e) 
no 



(18) 



where now e ~ Pujq- In the limit of low temperature, 
T ^ Wo, e becomes very large. It should be noticed that 
the factor In^t in Eq. (|l^) comes from the first double 
integral in Cuj{t), Eq. (25), whereas the factor Inwot in 
Eq. ( p^ comes from the next two single integrals of (25). 

12r 

10- 

o 



-1 



. 



1.5 



5 1 

Temperature T/ cOq 

FIG. 1. Temperature dependence of the factor e, given by 
Eq. (p^), in the limit of small friction, 7 ^ cjq. 

Discussion. From a technical point of view, the in- 
clusion of initial correlations between the system and the 
bath results in a modification of the integration contour 
in the complex-time plane (Keldish contour) that ap- 
pears in the path integral evaluation of the influence func- 
tional ]20|-|2^ . More precisely, the effect of the initial cor- 
relations is to couple the forward and backward integra- 
tion paths along the real time axis through an imaginary 
path along the Euclidean time axis, r = it. This leads to 
additional terms in the influence functional that are given 
by Euclidean integrals of the form J^^ (ir/(r). The pres- 
ence of the factor e in the propagating function (^ can 
be directly traced back to these terms. Clearly, even for 
very short times, when the dynamics of the system can be 
completely neglected, some of these terms, that depend 
only on temperature and not on time, will still be present. 
For initial thermal correlations, these terms explicitly de- 
pend on the nature of the system through the common 
initial equilibrium state of the system with the bath. As 
a consequence, even for arbitrarily small times, the atten- 
uation coefficient will be system dependent. Only in the 
special case of uncorrelated initial conditions is the atten- 
uation factor system independent at short times. It turns 
out, furthermore, that the contribution of these system- 
specific terms to the decoherence time will be negligibly 
small in the limit of high temperature, P 0, (as easily 
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seen from the Euclidean integral above) and will become 
very important in the opposite limit of low temperature, 
(3 ^ CO. This general discussion confirms the results ob- 
tained for the special examples of a free particle and a 
linear oscillator. 

Finally, it is also interesting to look at the temperature 
dependence of the factor e which can be written in the 
form 



47 



— r / duj 



LJC0th(/3cj/2) 



(^2 _ ujIY + 472^2 



(19) 



At high temperature, e asymptotically decays to zero as 
T^^, whereas at low temperature it diverges as T^^ [ggg 
Figure (|l|)]. In the absence of damping, Eq. ( p^ ) reduces 
to the simple expression 



coih{(3ujo/2)' 



7 = 



(20) 



We see from Eq. (|2y) that the high-temperature limit 
/9 — > is formally equivalent to the free particle limit 
Wo 0. This offers another explanation why the system 
dependent correction to the attenuation coefhcient be- 
comes vanishingly small at high temperature. Moreover, 
the presence of the hyperbolic cotangent in the denomi- 
nator of Eq. (^ ) hints at a possible connection between 
the divergence of e close to zero temperature and zero- 
point oscillations of the bath. 

In conclusion, we have examined the effect of initial 
correlations on the short-time decoherence of a super- 
position of two Gaussian wave packets. To this end, we 
have calculated the attenuation coefhcient for both a free 
particle and a linear oscillator, with and without initial 
thermal correlations. We have found that for factoriz- 
able initial conditions, the attenuation factor, and ac- 
cordingly the decoherence time, is system independent at 
very short times. On the other hand, for correlated ther- 
mal initial conditions, not only the temporal properties 
of decoherence are modified — changing from a stretched 
exponential to a Gaussian decay — but also the coher- 
ence length is affected. The latter is of the order of the de 
Broglie thermal wavelength. Moreover, the attenuation 



factor now has a system dependent term, containing the 
frequency of the oscillator and the temperature, and this 
even at times much smaller than the system characteris- 
tic time Ts ~ ^o"^- The system specific correction to the 
attenuation coefficient is small in the high temperature 
limit, where zero-point fluctuations are neglected, but 
becomes more and more important as the temperature is 
lowered. 

This work was funded in part by the ONR under con- 
tract N00014-01-1-0594. 



APPENDIX A 

In this appendix we collect for convenience the exact 
expressions of the quantities appearing in the propagator 
(^). The derivation of these quantities can be found in 
Rcf. The functions K(t) and N{t) are given by 



K{t) — V cotan vt 
^(,)_-exp[7t] 



sin vt 



(21) 
(22) 



where v — (tjg — 7^)^/2. We work in the limit of an un- 
derdamped oscillator where the friction coefhcient 7 is 
smaller than the frequency wq of the harmonic oscillator. 
For a free particle, — 0, one has to replace v i^. 
The variance of the position in equilibrium, k = ' 
is given by 



eg, 



27 



(a;2 



,2\2 



)2 + 472^2 



(23) 



Its asymptotic behavior at high (T^ujq) and low tem- 
perature {T <^ ujq) is respectively given by k ~ 1/i3ujq 
and K ^ Xj^ujQ. The functions C{t) and E[t) are, on the 
other hand, of the form 



dwtJC0th(/3w/2)/^(i) 



(24) 



where /a)(i) has respectively to be replaced by 



a(i) 



— \ — / ml dt"smiy{t~t')cosuj{t' -t") sin iy{t-t") exp[-f{t' + t")] 
sin vt Jo Jq 



E^{t) 



(C^2_^2)2. 


f 472^2 


sin vt 


2w(w2 _ 




1 


(C^2_^2)2_ 


f 4720)2 


sin vt 


2 


1 





•Smvt (Wq — (jj2)2 _|_ 4^2^2 



1 /■* 

/ dt! cos Lot! sin v{t — t!) exp[7t'] 

1 vt Jq 



{luI - uj-^Y + 472^2 ' 



dt' sin v{t — t') exp[7i'] ^(wg ^ cos uot' — 27a; sin wt'^ 



(25) 
(26) 
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APPENDIX B 



E{t) are of higher order and can therefore be neglected 
in the hmit considered in the paper]. We begin by com- 
In this appendix we evaluate the function C{t) = puting the first double integral in Eq. (25). Introducing 
sin^ vt C{t) [the contributions coming from the function the new variables u = {f + 1")/2 and v ^ t' - 1", we have 



a^{t) = / dt' dt" sin iy{t-t') sin iy{t-t") cos uj{t' -t")exp['y{t' + t")] 
Jo Jo 

du / dv sin iy{t — u — v/2) sini'{t ^ u + v/2) exp[2'-fu]cosuJV . 



10 JO 

The integral over u is readily obtained as, 



du sinv{t - u - v/2) sinv{t - u + v/2)e = 



2-yu 1 [ (e^''* — 1) COS i/w ^ 7Cos2t/i — i/sin2z^t — 76 
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7^ + i/^ 



(27) 



(28) 



(29) 



(30) 



at low temperature (the latter should be interpreted as a principal value p5[). Furthermore, by using the following 
two expressions. 



while the integral over v can be calculated using 

/ dcj (jj coth(/3cj/2) cosww ~ 2r j dcj cos ljv = 2ttT 6 (v) 

at high temperature, and 

/>oo 

du uj coth{f3uj/ 2) COS ujv / duj lu cos uv 



I dt' cosujt' sinv{t-t')e<^' 
Jo 

dt' sinujt' sinv{t ~ t')e^^ 



—v{ljq — Lo ) COS vt — ^{u + w ) sin vt + ve^^yiu^ — ur') cosujt + 27^ sin wi] 
2'^v cos vt^LD (7^ — + ) sin i/< + 75'''* [—27a; cos Lot + (wq — ) sin wt] 



we can rewrite the remaining three terms in Eq. (25) in the compact form 

2u: sin ut v sin vt exp[7t] — lJ^ sin^ vt 



The function C{t) is eventually given by 



(w^ - ^2)2 + 472^2 



C(t) = - / duJUJC0th((3uj/2)\aJt) + bjt)] 



(31) 
(32) 

(33) 

(34) 



[1] A.O. Barvinsky, A.Y. Kamenshchik, C. Kiefer, and I.V. 
Mishakov, Nucl. Phys. B 551, 374 (1999). 

[2] Ed. J. Wheeler, W. Zurek, Quantum theory and measure- 
ment, (Princeton University Press, 1983). 

[3] A. Galindo and M.A. Martm-Delgado, Rev. Mod. Phys. 
74, 347 (2002). 

[4] M. Brune et al, Phys. Rev. Lett. 77, 4887 (1996). 

[5] C.J. Myatt et al, Nature (London) 403, 269 (2000). 

[6] M. Mei and M. Weitz, Phys. Rev. Lett. 86, 559 (2000). 



[7] P. Mohanty, in Complexity from Microscopic to Macro- 
scopic Scales: Coherence and Large deviations, Edited by 
A.T. Skjeltorp and T. Vicsek (Kluwer, Dordrecht, 2001) 

[8] W.H. Zurek, Phys. Rev. D 24, 1516 (1981); 26, 1862 
(1982); Physics Today 44, 36 (1991). 

[9] D. Giulini, E. Joos, K. Kiefer, J. Kupsch, I.O. Sta- 
matescu, and H.D. Zeh, Decoherence and the Appearence 
of a Classical World in Quantum Theory, (Springer, 
Berlin, 1996). 

[10] G.W. Ford and R.F. O'ConneU, Phys. Lett. A 86, 87 
(2001). 

[11] G.W. Ford, J.T. Lewis, and R.F. O'ConneU, Phys. Rev. 

A 64, 032101 (2001). 
[12] D. Braun, F. Haake, and W.T. Strunz, Phys. Rev. Lett. 



6 



113] 



86, 2913 (2001). 
W.T. Stru nz, F. 
ph/0204129| ; W.T. 



ph/020510j 



Haake, and D. Braun, 
Strunz and F. Haake, 



^[uant 

quant- 



[14] V. Privman, |cond-mat/0203039| |quant-ph/0205037 

[15] J. P. Paz and W.H Zurek, in Les Houches, session 72. 
Coherent Atomic Waves, edited by R. Kaiser, C. West- 
brook, and F. David (Springer, Berlin, 1999), p. 533. 

[16] L.D. Romero and J. P. Paz, Phys. Rev. A 55, 4070 (1997). 

[17] R.P. Feynman and F.L. Vernon, Ann. Phys.(N.Y.) 24, 
118 (1963). 

[18] A.O. Caldeira and A.J. Leggett, Pliys. Rev. A 31, 1059 
(1985). 



[19] V. Hakim and V. Ambegaokar, Phys. Rev. A 32, 423 
(1985). 

[20] C. Morais Smith and A.O. Caldeira, Phys. Rev. A 

36,3509 (1987). 
[21] C. Morais Smith and A.O. Caldeira, Phys. Rev. A 41, 

3103 (1990). 

[22] H. Grabert, P. Schramm, and G.L. Ingold, Phys. Rep. 

168, 115 (1988). 
[23] G.W. Ford and R.F. O'Connell, Phys. Rev. D 64, 105020 

(2001). 



[24] F. Marquardt, cond-mat/0207692 



[25] A. Papoulis, The Fourier Integral and its Applications, 
(McGraw-Hill, New York, 1962). 



7 



